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SOLUTION OF SIMULTANEOUS LINEAR EQUATIONS AND COMPUTATION 
OF INVERSE MATRICES USING THE IBM CARD PROGRAMMED CALCULATOR 

Kenneth C. Rich 

U. S. Naval Ordnance Test Station 

China Lake, California 

Introduction 

The greatly increased use of matrix methods for solving certain problems en- 
countered in chemistry, aerodynamics, statistics and data analysis makes the devel- 
opment of machine methods for the solution of these problems highly desirable. The 
introduction of more complex computing machinery, such as the Card Programmed 
Calculator (CPC), which has certain definite characteristics and advantages over other IBM 
equipment, makes possible changes in procedure not previously incorporated in IBM 
methods. 

P. S. Dwyer reviews approximately 14 methods for the solution of simultaneous 
equations and inverse matrices by direct methods. 1 The method of single division as 
described by Dwyer appears to have better application to CPC operations than any of 
the others. Indirect or iterative methods at the present time may be performed on the 
CPC but do not appear to be as conveniently handled by it as by a direct method. In this 
paper a new direct method is presented which is an elimination method adapted from 
those suggested by Dwyer 1 and Frank Verzuh, 3 with modifications to fit the CPC. 

Form of the Matrix Equations 

A set of simultaneous equations may be expressed as 

anXi + a^x 2 + aisX 3 + . . . + amx,, = bi 
a 2 iX! + a 22 x 2 + a a aXa + . . . + aapc,, = b 2 



a^iXi + an2X 2 + a„3X3 + . . . + a^Xa = b n 

or expressed in matrix notation as 

AX = B 

1 Dwyer, P. S. "The Solution of Simultaneous Equations", PSYCHOMETRIKA, Vol. 6, 
No. 2 (April 1941), pp. 101-29. 

2 Verzuh, Frank. "The Solution of Simultaneous Equations With the Aid of the 602A Cal- 



culating Punch, " MATH TABLES, July 1949, pp. 453-62. 



By direct elimination methods we may solve for x^, and by back substitution solve for 
x»_ x , . . . , x 8 , Xi. Various existing IBM procedures use this scheme of elimination. 
However, if we introduce the composite matrix 2 * 3 



A :b 
-i : 6 



the elements of which are: 
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then by the elimination process we obtain a matrix of the form 
1 hi2 his . . . hi, h x 
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without the necessity of making a back substitution. The advantage of this method lies 
in the similarity of the arithmetical operation performed for all the computations. 

In a similar fashion we may compute the inverse by using in place of the (n + l)th 
column a unit matrix of n rows and n columns. In matrix notation this is expressed as 
1 B may be an n- columned matrix. For the examples given, however, B is a one- 
columned matrix. 



By direct elimination, the matrix 



a :i 
-i : o 



h : g 

• •••*• • 

o :c 



is now generated as follows: 
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where [C] is the inverse of [A]. This form appears to provide the quickest and most 
convenient method for CPC operation that we have found, proceeding as it does in a 
straightforward manner with identical computation procedures. 



Method of Reduction 

For the general case the reduction to the triangular matrix may be illustrated 
step by step by the following equations: 
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The advantage of this procedure is that rows 1A, 2A, 3A, 4A, . . . nA may be 
stored (up to the machine storage capacity), making it unnecessary to proceed with one 
elimination at a time, but practicable to operate as far as possible on one row before 
going on to the next. Card handling is thus reduced to a minimum, and elements of a 
particular row are handled only once. Since the rows of the matrix are read one at a 
time into the tabulator counters and since the tabulator counters accumulate, the only 
routine operations required by the CPC are those of multiplication, division, and transfer. 
Taking advantage of the form of matrix equations given above, the stored elements (1A, 
2A, 3A, . . . , nA) may be read back into the CPC, and by the identical routine operations 
the unknowns may be computed without the back substitution. 

By the above procedures, all row elements are read into the tabulator counters. 
When n > 6, operations are performed on the first portion of the row, after which the 
second portion of the row is read into the tabulator counters and the operation repeated. 
This is continued until all the elements of the row have been read into the counters. All 
the following rows are similarly handled. As the storage counters become loaded, inter- 
mediate results are punched out (in order) to be remerged into the following instruction 
deck. As an example of this procedure, the elements of a 10-by-10 inverse matrix are 
read into the counters first in a group of four elements, then in a group of six elements. 
In this way intermediate results are punched out twice (20 cards each) before the inverse 
is computed. The rows of a 20-by-20 inverse are split into a group of two elements, 
followed by three groups of six elements each. Nine distinct punching-out operations 
are necessary (80 cards each). The computing scheme will be illustrated in the following 
two examples. 



Five Simultaneous Equations 

The elements of rows 1-5 are read into the tabulator counters on instruction cards 
1, 7, 17, 30, and 45, respectively. Card Programmed Calculator operation codes will 
control the computation for each instruction card. All multiplications are read in 
negatively to quantities in the tabulator counters. (The instruction card numbers appear 
in the left-hand column. ) 

1 an a i2 ai3 a i4 a i5 bi 

o e a-12 ai 3 ai 4 ai5 bi > 

_ = K| =K| = KJ =K| __ --Kj K|. . . K| to storage 

a u au an 3-u an 

7 a 21 aaa aj 3 s.^ &as b2 

8-12 a|j al 3 a|4 ais b§ ala = a^ - K§a 2l ,etc. 

13-16 ais , a|* , ais ' bi , , 

—a- = L| —j- =LJ — = Lf —j- = L| Lf , . . Lf to storage 
a22 aaa a22 a 22 

17 asi a 32 ass a^ ass b3 

18-22 a| 2 a| s aJL a 2 x bf a§ 2 = a 32 - K|a 3 i,etc. 

23-26 ais alt a& b| a| 3 = als - L|a! a ,etc. 

27 " 29 _^=Mt _!* = Ml J?L = Mt Mi. . .M| to storage 

a3s ass ass 

30 a^ a*, a 4S a^ a^ b 4 . 

31-35 a.42 a 43 ali a«, b« aja = a 42 - Kla^ , etc. 

36-39 a|a al* ai b| a| 3 - a| 3 - Lia^etc. 

40-42 at a.^ b\ a^ = af 4 - M*af 3 ,etc. 

43-44 Jj = j^ 'Jd = N» N» , N| to storage 

a 44 a« 
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45 a 5 i as2 as3 as 4 ass bs 

46-50 a| 2 a| 3 af 4 ais b| ai 2 " = a 52 - K§a 5 i,etc. 

51-54 af 3 a| 4 af 5 b| af 3 = a§ 3 - L|a§ 2 ,etc. 

55-57 af 4 a£ 5 bt a.U = af 4 - M|af 3 ,etc. 

58-59 ais b| ais = af 5 - N|a| 4 , etc. 

= Ot Oe to storage 



60 



Js 



ass 



Note that down to this point we have triangularized the original matrix and have 
stored the elements of the resulting matrix. Here also we could begin the back sub- 
stitution, since 0§ is recognized as a solution. We prefer, however, to continue with 
the same procedure and operations, in effect reducing the triangular matrix to the dia- 
gonal form. 
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The entire procedure, from the time the first-row elements are read into the machine, 
takes approximately 45 seconds. 

Inversion of Fifth-Order Matrix 

The corresponding inverse, being, in a way, five sets of simultaneous equations, 
may be computed similarly, but at a considerable savings in time over five sets of five 
simultaneous equations. Instead of 102 instruction cards as needed by a set of simul- 
taneous equations, the inverse may be computed with 160 cards, in a time of approximately 
1 i minutes, by using the procedure indicated below in which twenty-five storage positions 
are required. (The instruction card numbers appear in the left-hand column. ) 
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To this point we have again triangularized the matrix, this time in the form of a diagonal 
band through the composite matrix. The inverse might be found by employing back 
substitution. Again, however, it is preferable to continue the routine operations as 
before and in this way to find the inverse. 
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The elements of the inverse, (1) to (5), are tabulated and summary-punched. 
Occasionally the diagonal element of a nonsingular matrix becomes very small, or zero. 
This condition is usually corrected by exchanging rows of the original matrix or by 
exchanging the rows after an intermediate summary-punching operation. Care should 
be taken not to change rows already being used in back substitution. These operations 
will have the effect of exchanging corresponding columns of the inverse matrix. If, by 
exchanging rows, the diagonal element remains zero for all combinations, the matrix is 
singular and no solution is possible. The presence of linear combinations of rows and 
columns is noted by a zero diagonal element. By inspection of the tabulation, the 
position of the corresponding rows and columns of the original matrix may usually be 
found. 

The examples given are indicative of the method with no restrictions as to the 
matrix being symmetrical or non-symmetrical. Expansion of the procedure for matrices 
of higher order has been done with great saving in time over previous IBM methods. For 



1 3 



example, by the method described in this paper, 4 a 10-by-10 inverse matrix may be 
computed in 11 minutes and a 20-by-20 inverse matrix may be computed in 90 minutes. 
Another IBM procedure using the CPC 5 requires only sixteen 10-digit storage units and, 
in addition, carries a check column. This procedure computes the 20-by-20 inverse 
matrix in two hours by use of Crout's method. This method has previously been judged 
the quickest IBM method available. The procedure described in this paper, however, 
has an advantage in that no additional equipment is required, such as reproducer, collator, 
or sorter. Hand merging of summary-punched cards behind marker cards has been 
found satisfactory in most cases and will not delay the computation time. Larger matrices 
require automatic merging of the cards. The use of forty-eight 10-digit storage units 
facilitates the reduction, but extra storage is not a necessary requirement. Changes in 
procedure make possible the same type of reduction using sixteen 10-digit storage units. 

The value of the determinant of the coefficients, if desired, is determined quite 
simply by obtaining the product of the leading or diagonal terms of each transformed set. 
These are expressed as 



A = An x Al 2 X AlaX A| 4 . . . A 



n 
nn 



Checks and Improvement of Accuracy 

The question of a check on the computations naturally arises. A check row or a 
check column may be carried along at all stages of computations. The sum of the new 
row elements or the sum of the column elements is carried as new columns or new rows, 
respectively. The operations performed on these are identical to the ones performed on 
the rows or columns of the matrix. 

4 A modification in the method of diagonalizing the triangular matrix has been suggested 
by B. G. Oldfield, one of the reviewers of this report. By making greater use of the 
Type 604 Electronic Calculating Punch for the accumulation of products, B - AC, the 
following saving in time may be made for the inverse : 

Number of cards saved = (n - l) 2 

where n is the order of the matrix. This saving in time is appreciable with large 
matrices. The cards saved by using the same procedure for a set of simultaneous 
equations will be 

n 

Number of cards saved = S (n - k) + (2n - 2) 

k= 1 

where n is the number of equations in the set. 

5 Chancellor, J. , and others. "The Solution of Simultaneous Linear Equations Using the 
IBM Card Programmed Electronic Calculator. " Proceedings, Industrial Computation 
Seminar (September 1950), International Business Machines Corporation, New York. 
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A check of the accuracy of final approximate results (Xo) may be made by computing 
the residuals 6 from the equation 

6 = b - AXo 
A second approximation, Xo + X x , may then be obtained by finding X x from 

AXi = 6 

Similarly, a check of the inverse as well as a more accurate inverse may be obtained by 
using the first approximation, Xo , to obtain Xo + Xi from the equation 

Xi = Xo (I - AXo) 

The latter two methods of checking and improving the accuracy require a large 
amount of computation which may be performed on the CPC but which is not conveniently 
handled. Therefore, checks of the computations are made by using the sum of rows or 
the sum of columns. The sum of the column elements carried as a new row will be in a 
better form for CPC computation. The addition of the check row increases the computation 
time very little. 

The method of solving simultaneous linear equations and the inversion of matrices 
described in this paper was set up specifically for use on the CPC. This procedure is 
believed to be the simplest and most convenient devised for IBM equipment. 
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6 Fox, L. , and others. "Notes on the Solution of Algebraic Linear Simultaneous 
Equations," QUART J MECH APPL MATH, Vol. 1 (June 1948), pp. 149-73. 
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NOTES ON THE CPC MODEL II AND COMPONENT MACHINES 

Eric V. Hankam 

Watson Scientific Computation Laboratory 
International Business Machines Corporation 



418 Pilot and Co-selectors 

1. I-p. u. of Pilot and Co-selectors 

a) If p. u. between 225* and 252°, selector will transfer, only for duration of p. u. 
pulse. 

b) If p.u. after 252°, selector will transfer immediately and hold until 225° of the 
next card, program or idle cycle. Any 418 cycle starts at 330° and therefore, 
if a selector is p.u. after 330°, it will hold until 225° of the same cycle. 

2. D-p.u. of Pilot Selectors 1-60 

Selectors hold through any intervening idle cycles in all cases below where selectors 
transfer until 225° (also applies to selectors 61-76). 

a) If p. u. between 225° and 235°, selector does not transfer at all. 

b) If p. u. between 235° and 252°, selector transfers immediately and holds until 
252°. 

c) If p. u. between 252° and 290°, selector transfers immediately and holds until 
225° (acts like I-p. u.) 

d) If p. u. between 290° and 350°, selector does not transfer at all. 

e) If p. u. after 350° and before 225°, selector transfers at 235° and holds until 
225° of the next card or program cycle. 

Note: Combining c) and e) above we see that if a p. u. pulse 
extends from before 290° to after 350° the selector will 
transfer immediately and hold until 225° of the next 
card or program cycle. 

3. D-p.u. of Pilot Selectors 61-76 

a) If p. u. between 225° and 235°, selector does not transfer at all. 

b) If p. u. between 235° and 290°, selector transfers immediately and holds until 
225° (acts like I-p. u.) 

c) If p. u. between 290° and 358°, selector does not transfer at all. 

d) If p. u. after 358° and before 225°, selector transfers at 235° and holds until 
225° of the next card cycle. 
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Note: By means of an externally pluggable "Hold" circuit for 
selectors 1-60, these selectors in case 2e above, can be 
dropped out like selectors 61-76 at 225° of the next card 
cycle, i.e. they will hold through any intervening program 
cycles. 

In the above description of selector timing, as well as in the notes to follow, an attempt 
has been made to supplement the timing charts contained in the Principles of Operation 
Manuals on the CPC and the 402 by including the exact timing of all such control panel 
hubs not described in these charts. 

4. 12-11 Pick Up Hubs 

The 12-11 p.u. hub does not accept on program cycles, but the 12-9 p.u. hub will 
accept impulses. 

5. Hold Circuits 

a) X-H hubs do not accept on program cycles. 

b) D-H hubs accept any impulse on card and program cycles. 

c) Statement b) implies that the Hold circuits of Pilot selectors 1-60 can be 
activated during programming in which case selectors will transfer for the 
remaining program cycles and the following card cycle (like selectors 61-76). 

d) The X-H and D-H hubs can be impulsed on any card cycle preceding program 
cycles. The Hold circuit will then be effective for any selector picked up 
during subsequent program cycles. 

Field Selector 

The Field Selector p.u. hub accepts only digit impulses; if impulsed by a "double" digit 
it wiU accept the first one, i.e. the larger digit. The selector transfers on the follow- 
ing card or program cycle. If Field Selector Hold is impulsed it will transfer on the 
following card cycle but, unlike Pilot Selector, it will not transfer throughout any inter- 
vening program cycles. 

Digit Selectors, Coding Selectors and Column Splits 

A coding selector can be transformed into an emitter on card and program cycles by 
wiring the DI impulse to the common hub of a coding selector. Coding selectors are 
operative on card and program cycles while Digit Selectors are operative on card cycles 
only. The DI impulse itself is active on card, program and idle cycles. Column splits 
are operative on card cycles only; at all other times the common hub is internally 
connected to the 0-9 hub. 

418 Counters 



Internal Connection of Counter Entries and Exits 

a) If a counter is impulsed to RI+, its entry and exit are internally connected. 

b) If a counter is impulsed to RI-, its entry and exit are internally connected. 
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c) If a counter is not impulsed to RI+ or RI-, its entry and exit are not 
internally connected. 

d) If a counter is impulsed both to RI+ and RI-, its entry and exit are not 
internally connected. 

2. Counter Read Out (RO) and Counter Read Out and Reset (RR) 

a) Impulsing counter TOTAL will cause counter RR. 

b) Impulsing counter TOTAL and RI- will cause counter RR. 

c) Impulsing counter TOTAL and RI+ will cause counter RO. 

d) Impulsing counter TOTAL and RI+ and RI- will cause counter RO. 

Combining 1. with 2. , we see that in cases 2a and 2d counter entries and exits are not 
internally connected while in cases 2b and 2c they are connected. Thus it is possible to 
interchange entries and exits in any problem. 

Whenever "Counter Exit Suppression" is impulsed, counter exits become inactive 
regardless of whether they are used as exits or entries. 

Statements 1. and 2. hold for card and program cycles. 

3. General Information 

a) If a counter is simultaneously impulsed to RI+ and RI-, digits will not 
read into counter entries. 

b) If an early All cycles pulse is wired to a counter entry position, the 
counter may not read in digits on the following cycle. 

c) Zero test: If the result of an accumulation in a counter is zero, all 
counter positions will contain 9's. Following a counter reset operation 
all counter positions will contain zeros. An X-pulse can be wired through 
an internal path between the CI and C hubs of a counter only when all 
counter positions are 9's. A RI- impulse to a counter following a reset 
operation will change zeros to 9's. 

d) Both counter exits and top counter exits emit an impulse for zero. 

e) Counter exits will read out the digit in each counter wheel and continue 
emitting all subsequent digits until zero is reached, e.g. a "3" will read 
out as 3, 2, 1, 0. 

f) Two or more counters should not be reset on the same cycle whenever 
their counter exits are split wired. 

g) If Counter Exit Suppression is impulsed by any digit (or 10-pulse) 
it will suppress the read out of all subsequent digits; e. g. if a 
6-pulse is used, digits 5, 4, 3, . . . will be suppressed. 
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418 Program Unit 

1. Special Program 

Whenever "Program Start" is impulsed, the feeding of cards is suspended and the 
machine goes through a series of program cycles until "Program Stop" is impulsed. 
The Channel entry, Minor, Intermediate, Major and Program 4 hubs form a 7 position, 
4-way "progressive" selector operative during special programming in the following 
manner: The seven Channel entry hubs are internally connected to the corresponding 
Minor hubs on the first program step, to the Intermediate hubs on the second program 
step, to the Major hubs on the third program step and to the Program4 hubs on the 
fourth program step. This internal connection is then repeated for the next four steps, 
etc. until "Program Stop" is impulsed. The internal connection lasts from 225° to 202° 
so that practically every impulse except 12-pulses can be wired through. 

The "program start" hubs (MI, IN or MA) should not be impulsed during Early All 
cycles time. The "program stop" hubs (MI, IN, MA) will accept from 320° on, but 
should not be impulsed later than 8-time. Programming can be started by impulsing 
any of the three (MI, IN or MA) program start hubs. Programming can then be stopped 
by impulsing a program stop hub of the corresponding or "higher" order. 

2. Program Couple 

Split wiring two or more adjacent "program couple" hubs will make corresponding 
program steps operative simultaneously. 

Split wiring two non-adjacent "program couple" hubs will make the two corresponding 
program steps operative simultaneously, and will skip the program steps in between. 

(Note the similarity with the 602 -A program unit) 

If program steps are repeated by means of the special program unit, the four pro- 
gram couple hubs cannot be splitwired arbitrarily. They can be "coupled" in the 
following combinations: 



The remaining combinations fail: 



1, 


2 




2, 


3 




1, 


2, 
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1, 


4 




2, 
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3, 
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1, 


2, 
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Run-in and Run-out Operations 

1. Last Card Auto Total Switch 

a) This switch, if on, causes a Major program start on the run-in and run- 
out of cards. Furthermore, it suspends the emission of impulses from the 
comparing exits for any purpose, but it will not suspend programming 
unless program start is impulsed from comparing exits. 

b) Whenever the switch is on, and the special program feature is used, 
programming must be stopped by impulsing the "Major Program Stop" 
because the machine would otherwise never stop programming on the run- 
in when a major program step is started by the switch. 

2. LC Hub 



a) The LC hub emits a pulse (162° to 217°) on the run-in, i.e. when the first 
card is at the 2nd Reading Station, and on the run-out, i.e. when the last 
card is at the 3rd Reading station. 

b) If the Last Card Auto Total switch is on, the LC hub can be used as inlet 
from 320° to 290° on any card cycle. Its function is to cause a Major Pro- 
gram start. 

3. Final Total 

A counter will not convert on the run-out when its TOTAL hub is wired from Final 
Total. However, it can be made to convert if the TOTAL hub is wired from an All cycle 
instead of the Final Total through a selector picked up on the last card. 

Final Total conversion will be added as standard feature on forthcoming 41 8' s and 
can be installed on earlier models. 

4. Printing 

Printing and spacing is automatically suspended on the run-in, i.e. on all cycles 
prior to the card cycle of the first card at the 3rd Reading brushes. 

Hammerlocks 



The Hammerlock I-p.u. can be used on card and program cycles, but the D-p.u. , 
if impulsed on a program cycle, will operate on the first card cycle after programming. 

Stop Controls 

There are two hubs on the 418 panel which, if impulsed, will cause the machine to 
stop. The "Test" hub accepts any impulse, while the "Stop" hub accepts only during 
Early All cycles time. 

In CPC operations the "Stop" hub will stop the machine after the bus has completed 
its calculation. If the "Test" hub is impulsed not later than 0-time the 605 will not 
start programming on that same 418 cycle, and cannot be single cycled. If the "Test" 
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hub is impulsed at X-time the 605 will not start programming but can be single cycled. 
(This applies to later model machines where 605 calculate time has been advanced 5°, 
i. e. where calculate time extends from 195° to 280° instead of 190° to 275°. ) In both 
cases the 605 counter will not reset if wired to do so on that 418 cycle, i.e. the counter 
will still contain data from the preceding calculation. If the "Test" hub is impulsed at 
12-time the 605 will complete its calculation before it stops. 

Summary Punching 

The summary punch PU hub accepts any pulse except in the range from 330° to 352°. 
A summary punch (idle) cycle will occur at 330° following the impulse to the PU hub. 
Therefore, if the PU hub is impulsed by an All cycles, summary punching will take place 
right after Early All cycles time; on the other hand if a 10- or 9- pulse is used to 
activate the PU hub, summary punching will not take place until the end of that cycle at 
330°. * 

10 Pulse Timing 

On some 418's the timing of pulses occurring at 10-time is different from that 
indicated on the CPC timing chart. The 10-common emits from 349° to 3° while the 
10-hub, the * symbol and the 10 segment of the credit symbol exit emit from 348° to 2°. 

Channel C Shift 

The 15 position Channel C Shift Entry and the 10-position Channel C form a "shift" 
selector that can be used on the 418 independent of the CPC; if desired. For all 
operations, it is controlled by the "Shift" (pick-up) hub which will accept 1,2, 3, 4 and 
5 pulses, as well as multiple pulses 1, 5 or 2, 4 or 1, 2, 3 (each adding up to 6). Let 
us denote the 15 position of Channel C Shift entry by 15, 14, ... 1 from left to right. 
Normally, positions 15 to 6 are internally connected to the Channel C hubs. Channel C 
is internally connected to positions 14 to 5 if the "Shift" is picked up by a 1, to positions 
13 to 4 if the "Shift" is picked up by a 2, etc. up to a shift of six positions. This 
internal connection lasts from 235° (following the "Shift" pick up) to 202° of the following 
card or program cycle. 

Associated with the "Shift" hub on the 418 control panel and also under its control, 
is a 1-position 6-way Field selector (called Channel shift) on the 605 control panel. 
Normally its internal connection is from C to O; if the "Shift" is picked up by a 1, 2, 
... 5 the selector will transfer to the correspondingly labeled levels during 605 
calculate time. Similar to a Calculator Selector, this field selector can be used for any 
selective operation such as selection of program pick-up or suppression. 

Channel A (or B) Card Read In 

A zero wired to the left position of the Channel A (or B) control hubs will establish 
an internal connection between Channel A (or B) and Channel A (or B) Card Read In. 

* NOTE: Digit and X pulses should not be wired to the summary punch PU hub if 
summary punching is to take place on two successive 418 cycles; 12 pulses can be used. 
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The adjacent sign hub (not to be confused with the "Channel Minus" A (or B)) provides 
an eleventh position of this "selector". The transfer time is from 235° following the 
pick up by the zero pulse to 202° of the following card or program cycle unless "Card 
Feed Channel Hold" is impulsed in which case it will transfer on the first card cycle 
following programming. 

Counter Channel Control 

The right (units) position of the Channel A, B and C Control hubs determines which 
of the 8 rows of the Counter Channel Control hubs become active. Multiple digit pulses 
cannot be used to activate more than one row at a time. Only the largest digit will be 
recognized. 

605 Program Pick up 

a) Each of the three program exits of a single program step emit isolated pulses. 
However, the same source pulse is operative for corresponding exits of series 
A, B and C, i.e. corresponding exits of series A, B and C emit the same 
impulse. Therefore, backcircuits will usually result if the same block of 3 
program steps is picked up simultaneously for more than one series. 

b) It is not possible to pick up the same block of 3 steps for more than one series 
even though mutually common steps are suppressed so that no step in that block 
is active on more than one series at once. The reason for this is that the 
suppression pulse makes the internal source pulse which supplies the program 
exit pulses to series A, B and C inoperative. 

Program Suppression 

a) Do not split wire two different suppression pulses into the same program 
suppression entry hub; always use selectors or filters. 

b) Suppression pulses should be wired only to program suppression entry hubs. 

For example, a suppress on zero pulse should not be wired to a group suppression 
pick up. (To accomplish this, wire a program exit, which is suppressed on non 
zero, to the group suppression pick up. ) 

605 Zero Test and Negative Balance (N. B. ) Test 

1. N. B. Test 

A N. B. test does not change the number tested in the electronic counter (EC). 
sSuch a test can be made on any program step except during multiplication or 
division operations. In particular, it is possible to impulse EC HI, N. B. test 
and EC RE on the same step. 

2. Zero Test 

a) A Zero test will, in general, change the number tested in the electronic 
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counter. No other operation should be performed in the counter together with 
Zero test. The counter can be reset on the same step, however, since a 
counter reset operation does not occur until the beginning of the next program 
step. The shift unit can be used, and numbers to the right of the position tested 
will be ignored. If the number tested is positive, an internal source will enter 
a - 1 into the position tested. If the number tested is negative a •> 1 will first 
be entered into the units position without carry, and then a + 1 will be entered 
into the position tested with carry. Thus, because of the loss of carry, negative 
numbers ending with a 9 in the units position will be altered differently from 
other negative numbers following a zero test. 

All zeros are treated as positive numbers with the exception of a zero remainder 
following a division with a negative dividend. (If desired, it is possible to 

change such a - remainder to + in the counter simply by impulsing EC RI+ . ) 

* 

b) As seen above, an internal source emits a 1 into the position tested with the 
opposite sign of the number in the counter. If the tube (located at panel 1-5 U) 
that emits the 1 were removed, any digit besides a 1 could be emitted from the 
emitter on the control panel into the counter on a zero test. Also, if a number 
were transferred from a storage unit to the counter and zero test is impulsed, 
it would go into the counter with the opposite sign of the counter (regardless of 
its own sign). If the tube (located at panel 3-5H) that causes the units position 
of a negative number to be decreased by 1 (it is not an actual 1 pulse) were 
removed, any number except -1 could still be zero tested; but a - 1 in the 
units position would be recognized erroneously as a zero. 

In the 1/2 Adjust circuit a tube emits a 5 into the counter with the sign of the 
number standing in the counter. If this tube (located at panel 1-7T) were re- 
moved, any number, regardless of its sign, could be read from a storage unit 
into the counter under the sign of the counter by impulsing 1/2 Adjust. This 
would greatly facilitate manipulations with absolute values. For a normal 1/2 
Adjust operation it would be necessary to wire a 5 from the emitter. 

c) It is not possible to Zero test on a program step which is suppressed on zero as 
a result of a previous zero test. 

605 Selectors 

1. Calculator Selectors 



a) On CPC operations, calculator selectors can be picked up (on the 418 
control panel) from 350° to 310°. The selector transfers immediately 
until 310°. It is essential, however, that the pick up pulse be effective 
before 605 calculate time (190° - 275°). Thus, a 418 X-pulse will be too 
late to transfer the selector before the 605 starts programming. Incidentally, 

* The tube removals mentioned in this section can be handled by special arrangement 
with the IBM Corporation. 
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such a condition cannot be detected by single cycling the machine since in 
that case the selector will transfer before the first program step is single 
cycled on the 605. 

b) A Program source pulse may be wired through the transferred points of a 
Program Repeat Selector to the coupling exit of a Calculator selector. In 
this case the coupling exit serves as entry to pick up the selector (it is not 
common with the pick up hub on the 418 control panel). The calculator 
selector will then pick up at the same time as the Repeat Selector and 
remains transferred until the end of the 605 calculate time unlike the 
Repeat Selector which may be dropped out at the end of any program sweep. 

c) A coupling exit of a calculator-selector p. u. on the 418 control panel can 
be wired to the coupling exit of another calculator selector which will then 
transfer together with the first selector. 

d) The pick up hub of a calculator selector on the 527 control panel will not 
accept a pulse before X-time, i.e. a Y-pulse or a 527 Pilot Selector 
coupling exit pulse cannot be used directly to pick up the selector. (This 
differs from the 521-604 machines). 

Program Repeat Selectors and Repeat Delay 

a) A Repeat selector, once picked up, will ordinarily drop out at the end of 
the earliest program sweep when Repeat Delay is operative (or the end of 
605 Calculate time, at the latest). However, if the selector is picked up 
again it will remain transferred (for at least another program sweep) 
regardless of whether Repeat Delay is operative or not. 

b) If Repeat Delay is both picked up and dropped out on the same program 
sweep, the later instruction will be effective. If it is picked up and dropped 
out on the same program step, the outcome is indeterminate. 

c) Storage assignment can be selected through a Repeat Selector. Storage 
assignment exits and Program source pulses are internally common, they 
emit pulses throughout 605 calculate time, but are broken at the end of 
each program sweep whenever Program Delay is operative. Calculator 
coupling exit impulses, on the other hand, are continuous pulses throughout 
calculate time and should not be wired through a Repeat Selector. (Being 
continuous pulses, they would arc selector points.) They can, however, 

be wired directly to program pick up or storage assignment entry hubs. 

d) Method to obtain 180 distinct program steps: 
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Prog. p. u. Series A ■< o N 

Prog. p.u. Series B^ o T 



Repeat 
Selector 1 



N 

Prog. p. u. Series C •< o T Re P eat 

Selector 2 
C 



Prog. Source- 



Program Repeat, Repeat Sel. 1 
p.u., and Repeat Delay p. u. are 
impulsed by series A program 
exits. (Any of the 60 steps can 
be used and all £hree instructions 
can be wired from the same step 
if desired. ) 

Program Repeat, and Repeat sel. 
2 p. u. are impulsed by series B 
program exits. 



Thus, the machine will go through three sweeps A, B and C in that order. 
If program repeat is impulsed again on series C, the machine will repeat 
the 180 steps until any of the Program Repeat steps is suppressed. It is 
possible to obtain more than 180 distinct steps by picking up program 
repeat selectors (to which Calculator selectors can be "coupled" as 
described in the section under Calculator selectors) at the end of sweep C, 
and by wiring program exits of series A through these selectors. 

In another variation of the problem, if Repeat Delay d. o. were wired on 
series C, Repeat selector 2 would remain transferred, and the machine 
would now start with series A, B, C and then repeat C until the Program 
Repeat instruction on that series is suppressed. 

605 Channel and Shift Unit 



a) Whenever two distinct digits are transmitted along the same channel position on 
the same step (e.g. if two storage units are simultaneously impulsed to read 
out into the counter), the result will be the absolute value of their difference 
without carries, and a minus sign if either of the two original numbers were 
minus. 

b) It is possible to emit a single digit from the Emitter into the 2nd, 3rd, . . . 6th 
position simultaneously by impulsing the corresponding shift unit on the same 
program step. 

CPC Sign Control Timing 

There exist two equivalent representations for negative numbers in CPC calculations. 
Negative numbers read out without conversion from 418 counters into the 941 will appear 
in complement form recognizable by a 9 in the eleventh position of a 941 storage unit. 
Negative numbers read out from the 605 into the 941 will enter the 941 as "true negatives", 
i. e. , as a true figure with a 5 in the eleventh position of a 941 storage unit. 
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1. 418 Counter N. B. Test Exit 

The N. B. test exit emits an impulse (258° - 270°) on card and program cycles, 
whenever the counter contains a negative number and the NB - AC switch is wired on. 
(If the switch is not wired, the NB hub emits only on program cycles. ) A NB test exit 
of a counter wired to the Negative Balance Control of the same or another counter will 
cause a conversion (idle) cycle in the latter counter whenever its TOTAL hub is impulsed. 
Except for printing purposes, it is almost always desirable to eliminate 418 conversion 
cycles during CPC operations. For example, it will save time to read from a 418 counter 
without conversion to the 605. A complement read out of a 418 counter can be converted 
to a true figure with a sign in the process of reading it into the 605, i.e. , a 605 storage 
unit will register the complement of each digit read into it whenever the sign hub of that 
storage unit is impulsed between 235° and 275°. The N. B. test exit pulse falls in that 
region. 

2. Top Counter Sign Hub 

This hub emits an impulse identical to the N. B. test exit at 258° - 270° on card and 
program cycles (regardless of whether the NB - AC switch is on or off). 

3. 418 Credit Symbol Exit 

The Credit Symbol exit of a counter emits an impulse on a card or program cycle 
following a conversion cycle of that counter. The NB - AC switch must be on for this 
hub to emit on card cycles. It also emits an impulse when a number is read into the 
counter negatively on a card cycle. The impulse consists of two parts: 282° to 296° and 
a 10 pulse. 

4. Channel A (or B) - (control panel location BG-BI, 11-12, 23-24) 

a) If a complement is transmitted along Channel A (or B) from the 941, the " - " 
hub emits from 258° - 270° and a 9 pulse. 

b) If a true negative is transmitted along Channel A (or B) from the 941, the " - " 
hub emits from 282° - 297° and a 10 pulse. 

c) If the Channel Minus A (or B) hub (control panel location AO: 55-58) has been 
impulsed by an X- or Y- pulse, the " - " hub will subsequently emit from 282° 
- 297° and a 10 pulse. 

d) If the " - " hub is impulsed by a NB test exit (or Top Counter sign) it will emit 
a 9 pulse, i.e. , in this case the " - " hub is both an inlet and outlet. 

It should be noted that the main function of the 258° pulse (if wired to a 605 storage 
sign hub) is to transform complements to true figures upon entry into 605 storage units. 
The 282° pulse can be used to activate 418 counter sign reversal. The 10 and 9 pulses 
are useful in printing operations. 
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5. Channel C - 

If a number is transmitted to a 941 storage unit via Channel C, and 

a) the " - " hub is impulsed by an X- or Y- pulse on the previous cycle, or by a 
10<- pulse, a "5" will register in the eleventh position of the 941 storage unit, 
identifying a true negative. 

b) the " - " hub is impulsed in the region from 235° - 270°, a "9" will register 
in the eleventh position of the 941 storage unit, identifying a complement. 

The " - "hub will not accept a 282° to 297° pulse. 

Digit pulses to 9 must not be wired to the " - " hub. 

6. FS, GS, and MQ Sign Entry 

Accepts any pulse from 255° on to store a minus sign. The sign trigger itself does 
not turn on to register the sign until 182°, providing the corresponding storage unit is 
also impulsed to read in. The sign hub should always be impulsed before 605 calculate 
time, i. e. , with an 0- pulse at the latest.* The sign hub has a dual function. As previously 
explained, if impulsed between 235° and 275°, it controls the subsequent read in of digits 
so that the storage entry will register the complement of the digits. 

7. GS and EC Sign Exit 

Emits at 305° - 320° and a 10 pulse whenever a negative number is read out of the 
corresponding units. 

8. Counter Sign Reversal 

Accepts any pulse from 270° to 327°. 

9. Channel 9, A (or B) 

Any impulse from 235° to 270° wired into Channel A (or B) -, or developed in 
Channel A (or B) - hub (e.g., when transmitting a complement from the 941) will make 
Channel 9 hubs emit 9 pulses. 

Zero Test 

The zero test hub on the 418 Control Panel will emit a 10 pulse following a 605 
calculation if during that calculation any Zero Test indicated a non-zero condition. 

Transfer Operations 

a) Whenever the 418 Counter Reconversion feature is used, a counter should not be 
read out on two successive card or program cycles, since the reconversion 
cycle would coincide with the conversion cycle of the second read out. The 
second read out would then be in complement form. 

b) A number should never be read out from the 941 on the cycle after it has been 
entered into the 941. At least one 418 cycle should intervene between entry and 

*If the sign hub is impulsed earlier than 330° and an idle cycle occurs, the sign will not 
hold through that idle cycle. Thus, to preserve the sign, the sign hub should be impulsed 
after 330°. 
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read out. In "standard" CPC operation this will mean that two instruction 
cards must pass through the machine between read in and read out since instruc- 
tions that control the reading into the 941 are frequently wired from the 3rd 
reading station while read out instructions are wired from the 2nd reading 
station. 

c) The 941 can be used as auxiliary storage for the 418 independent of the 605. The 
ST - INT switch does not have to be wired on in order to read into the 941. It 

is merely a safeguard since the machine will not run if the switch is wired and 
the 941 power switch has not been turned on. 

d) It is possible to read one number via Channel C into more than one storage unit 
of the same 941 bank. This is accomplished by impulsing the units position of 
Channel C Control with any combination of the digits 1, 2, . . . , 8 which identify 
the corresponding storage units of a 941 bank. 

e) If the 605 is used by itself, numbers read from 527 brushes cannot be entered 
into a 605 storage unit unless the storage unit is impulsed to read in. On the 
other hand, in CPC operations, the entry of numbers into the 605 is wired on 
the 418 control panel and numbers can enter 605 storage units without a read in 
impulse. However, unless the storage unit is impulsed to read in, it will not 
clear, and numbers will enter cumulatively without carries. 

f) If two 605 storage units, say Gl and G3, are splitwired to RO on the 605 control 
panel, and only one of them, say G3, is impulsed to RO on the following 418 
cycle (wired on the 418 control panel), Gl would also receive an RO impulse. 
This would cause an error if Gl were wired to RI on the 418 control panel. 

g) If a negative number is read out of the electronic counter (during a 418 cycle) 
and EC RR is wired from a 418 All cycles (or Card cycles) instead of an 
Electronic All cycles (or Electronic Card cycles), the number read out will be 
less by 1 in all positions. 

h) Do not split wire Channel A (or B) common or 10-common to any source pulse. 
This may cause the erroneous entry of a number into the 941 through an internal 
back circuit. 

527 Features During CPC Operations 

a) The following hubs on the 418 control panel are common with correspondingly 
labeled hubs on the 527 control panel: 

Calculate selectors pick-up hubs 

FS, GS, and MQ entries 

FS, GS, and MQ sign entries (not common on idle cycles) 

FS, GS, and MQ RI are common at 10 time only 

GS and EC exits 

GS and EC sign exits (not common on idle cycles) 

* Also, a storage unit in the 941 should not be impulsed to read in on two successive 
418 cycles. 
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GS RO and EC RO and RR 

The twelve SP Control Entry hubs on the 418 control panel are common 
with the 12-11 hubs of the 12 column splits on the 527 control panel. 

Any of these hubs (providing they are not used for their regular function) can be 
used as "connectors" to transmit an impulse between the 418 and 527. For example, a 
527 Pilot or Punch Selector can be picked up under the control of the 418 by wiring a hot 
527 pulse to the 418, through a 41& latch selector, back to the 527 and the Pilot or Punch 
selector pick up hub. If the hot 527 pulse is a Y-pulse, and if it is wired to the immediate 
pick up hub of the pilot or punch selector, summary punching of digits and X-pulses can 
be selected on the same cycle. 

b) Calculator selectors cannot be picked up with 527 pulses. Negative Balance 
Selectors and Zero Check on the 527 control panel cannot be used in CPC 
operations. The DPBC and Overflow features are operative in CPC operations. 

CPC Channels 

The block diagram illustrates the principal internal channels on the CPC. Except 
for summary punching, all pluggable exits and entries appear on the 418 control panel. 
Any exit can be wired to any entry to complete a path for transmitting numbers between 
the various units. For example, it can be seen from the diagram that a number from the 
941 via Channel A can be wired to any of the 605 storage entries, 418 counters, type 
bars or even via channel C to another 941 storage unit; but it cannot be wired to the 527 
punch magnets directly. Another possibility, evident from the diagram, is the connection 
of a 605 exit to another 605 entry. Computations on the CPC frequently require the 
saving of a result obtained in the 605 for the following 605 calculate time. This can 
easily be accomplished on the 605 control panel, but the idea suggested above might save 
badly needed 605 program steps. 

The author wishes to thank Mr. E. Lupano, IBM Customer Engineer, for his 
assistance in the preparation of this article. 
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ACCOUNTING MACHINES, TYPES 402 - 403 



TIMING CHART 
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A CARD-NUMBERING SYSTEM FOR CPC CALCULATION 

William C. McGee 

General Electric Company 

Richland, Washington 



Have you ever stood at the CPC wondering which deck of cards to feed in next? 
Have you ever thought to yourself, "Now, did I or didn't I change those loading cards?". 
Have you ever dropped a CPC instruction deck and had to start again from scratch? If 
you have, you would probably welcome some scheme for avoiding these situations which 
have a tendency to repeat themselves with annoying frequency. 

With just such situations in mind, we have developed a card-numbering system 
for CPC calculation which has successfully eliminated almost all of the card-handling 
problems at the CPC. 

Aside from identifying the cards and keeping them in order, you may ask why a 
card-numbering system is necessary at aU. After all, if your calculation requires 
1000 cards, why not number them consecutively from 1 to 1000? In this manner they 
will be positively identified and should you drop them, very easy to sort back into order. 

The answer to this question is found in the concept of the loop . A loop is simply 
a part of a calculation which is repeated a certain number of times. Each time it is 
repeated, new data is supplied to give new results, unless redundancy is desired, e. g. , 
in checking. On stored-program calculators, executing loops (i.e., getting into them, 
performing them, and getting out of them) is automatically taken care of by the 
instructions loaded in the memory. On the CPC, however, loops must be executed by 
means of card handling alone. This means that if a certain part of the calculation is to 
be repeated, say 5 times, the deck of cards performing that part of the calculation must 
be fed through the machine 5 times, each time with different data supplied. 

Loops may vary in length. Loops which are very long (say 1000-2000 cards) are 
probably most efficiently executed by hand, i. e. , by running the whole loop through as 
a unit, each time manually changing the data cards. This is not efficient for shorter 
loops, especially if they are to be repeated a large number of times, because of the 
excessive card handling it requires. In the past, shorter loops were executed by 
reproducing the cards in the loops as many times as the loops were to be repeated, and 
swapping in the data cards, either manually or on the collator. This was the right 
approach, but suffered from the fact that the new cards created by reproducing were not 
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uniquely identified. If the deck were dropped, or part of it returned to the wrong place 
in a tray, the jig, for all practical purposes, was up. 

The main purpose of this card -numbering system, then, is to identify uniquely all 
cards made for the purpose of executing loops in CPC calculation. It is interesting to 
note that if the problem involves no loops, or the loops are to be executed manually, then 
no card-numbering problem is involved. The cards may simply be numbered sequentially. 
If this situation obtains, the system to be described below reduces, as it should, to one 
of simple sequential numbering. 

The Flow Chart 

The first step in preparing a deck of instruction cards for CPC calculation under 
this system is to draw up a flow chart of the computation. A flow chart is simply a 
summary of the sequence of operations in a given computation. It serves the same 
purpose as a formal outline of computational procedure, but has the advantage of being 
concise and graphic. It is valuable for two reasons: (1) it locates each operation in the 
calculation, and (2) it defines each loop to be used in the calculation. Both of these 
features will be indispensable when we come to numbering the cards. 

The flow chart is constructed from three distinct elements which we have called 
(1) Operations, (2) Decisions, and (3) Modifications. 

Operations: An Operation is defined as anything that a card can do. Whether it is 
an arithmetic operation, obtaining a function, making a branch test (if such is available), 
modifying a variable by adding increments, punching, printing, or restoring the paper, 
it is called an operation if it can be accomplished by means of codes punched in a card. 
It may appear that this definition is somewhat broader than that usually given for a flow 
chart operation, in that it includes the branch test (a decision element) and a special 
type of modification (i.e. , by adding increments). It turns out that unless we consider 
these special card operations as Operation elements on the flow chart, confusion will 
arise over their interpretation. 

Decisions and Modifications : In contrast to the Operations, it might be said that 
Decisions and Modifications are anything that a card cannot do. Decisions must 
generally be made at the end of the loop; the decision will be either to go back through 
the loop again, or to proceed with the rest of the calculation. This decision will always 
be of the unconditional type, i. e. , not depending on any calculated result. (Recall that 
conditional decisions on the CPC can be made only with the branch test device which we 
have agreed to consider an Operation. ) Modifications are usually made after a decision 
for the purpose of making the next traverse of a loop unique. Again, this type of 
modification is the "card swapping" type, and should not be confused with the "incre- 
mental" type described under Operations above. 
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It is important to note that the decisions and modifications as defined here are 
implied in the manner in which the instruction deck is assembled. The programmer 
does all the deciding and modifying necessary to the calculation when he puts the deck 
together. The codes punched in the card have nothing to do with Decisions and 
Modifications. 

With the elements of the flow chart defined, we are now in position to draw a flow 
chart for any given problem. Let us take for an example a problem in numerical inte- 
gration. Suppose we wish to prepare a table of Y(a, b t ) versus b 4 where Y is defined 
as follows: 

Y(a,bj) = / FfcbOdx, 1=1,2,3,4,5. 

o 

Using Simpson's rule, we may write approximately 

Y(a, bO = ^- { F(x„ , bO + 4[F( Xl , b.) + F(x 3 , bi) + . . . + F(x 9 , b,)] 

+ 2[F(x 2 , bi) + F(X4 , bO + . . . + F(x 8 , b t )] + F(x l0 , b s )} 

where we have divided the range of integration to a into ten intervals. The flow 
chart for this computation is shown on Chart 1, To emphasize the distinction between 
functions performed by cards (the Operations) and by card-handling (Decisions and 
Modifications), we have placed all the Operations in a vertical column, and all the 
Decisions and Modifications to the side of this column. 

Programming 

The next step is to number the Operations consecutively, starting from the top. 
In our example, there are nine Operations. The numbered Operations are called blocks, 
and the numbers assigned to them are called block numbers . Each block represents a 
single independent part of the calculation, and may thus be programmed as a single 
independent unit. The programming for a block may require any number of cards, even 
just a single card. Each card in a block is given a sequence number to establish its 
position within the block. Blocks consisting of one card only are given the sequence 
number zero. 

At this point each card programmed has a block number and a sequence number. 
On our CPC systems we have reserved one column (column 1) for the block number and 
three columns (columns 4-6) for the sequence number. This implies that we will never 
encounter problems with more than nine blocks, or more than 999 cards in a block. The 
allotment of columns for these numbers is, of course, entirely arbitrary. We could 
just as well have reserved two columns for the block number and two columns for the 
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sequence number, giving us a maximum of 99 blocks, with 99 cards in each block. In 
fact, it is quite possible to leave this matter flexible, to be determined by the needs of 
each problem. The subsequent procedure will, of course, have to be modified slightly 
for each problem. 

A third number is required now to identify the cards in the programming process. 
This number is called the set number and is used to distinguish cards which perform 
the same operation (e. g. , loading or spread entering) on different sets of data. The 
set number is really synonymous with the "modification number" and its need is clearly 
signalled by the presence of Modification elements on the flow chart. The modification 
is understood to apply to a whole block, rather than a single card, even though many 
blocks requiring modification may consist of only one card. The blocks requiring 
modification, and therefore requiring set numbers, are those immediately following a 
Modification on the flow chart. Referring to our example, it is seen that Block 1, which 
is a load card, follows the modification element in the outermost loop and therefore 
requires set numbers. Since bi will be modified five times we assign five set numbers 
to Block 1. It is convenient to refer to blocks requiring modification by their Block-set 
number, and in our example we must prepare cards for Block-sets 101, 102, 103, 104, 
and 105. In our CPC system, two columns (columns 2-3) are reserved for set numbers, 
allowing for a possible 99 modifications of a given block. This again, is arbitrary. A 
majority of the blocks will not require modification, and these are given the set number 
zero. 

Note that there is another modification in our example which we appear to have 
ignored, namely, the stepping of x by the increment (a/10). Had we decided to step 
x by loading a different value of x each time, we certainly would have to set-number 
Block 4. But we have chosen to do this as a card operation, and purposefully omitted 
the modification element so that no additional set-numbering is called for. 

The cards in the instruction deck programming now bear three numbers, arranged 
(arbitrarily) as follows: 



BLOCK 
NO. 


SET 
NO. 


SEQUENCE 
NO. 


1 


2 3 


4 5 6 



A sample of the programming is shown on Chart 2. Note that the entries on the pro- 
gramming sheet are not necessarily in the order that the cards they represent will pass 
through the machine in solving the problem. Rather, they are in Block- set order. This 
makes it possible to group jd^ of the information that will be used in the problem, which 
is extremely useful in checking back to the original data, and somewhat more 
aesthetically appealing than carrying this data around on the back of a card or in one's 
head. 
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Checking the Program 

It is usually desirable to test the programming after it has been keypunched. Since 
the deck has not been expanded, this process will require careful attention from the 
person doing the testing. However, since the test is ordinarily not too long, the con- 
centration span of the tester is not too protracted, and cards fed in the wrong order will 
not incur a serious loss. The lack of a procedure for testing the program is probably 
the weakest point in this system. 

Expanding the Instruction Deck 

When the program has been checked out the instruction deck is ready to be expanded 
for the production run. At this point it is necessary to introduce a fourth and final card 
number which we call the group number. Group numbers are required to sort the 
expanded deck into the order in which it will pass through the machine. They are deter- 
mined with the aid of the flow chart as follows. 

The first group begins at a (symbol denoting start of computation), and ends when 
the first modification is encountered (either a card modification or a card-handling 
modification). The second group begins where the path from this modification intercepts 
the operation column and ends when a modification is again encountered. In this manner 
it is possible to number all groups, up to and including the last group, which ends at /3 
(symbol denoting the end of computation). 

For example, in our problem the first group begins at a , and ends at the first 
modification, which happens at Block-set 400 after the first time around the inner loop. 
Thus Group 1 consists of Block-sets 100, 200, 300, 400, 500, and 600. Group 2 begins 
with Block-set 400, and ends after the second time around the inner loop. Thus Group 2 
consists of Block-sets 400, 500, and 700. This process continues up to Group 10, which 
consists of Block-sets 400, 500, 800, and 900. This information is summarized on the 
right-hand side of the CPC Deck Assembly Request (Chart 3). 

Strictly speaking, our group numbering is not complete, since we have not reached 
/3 . Had we done this, we would have 50 groups instead of 10 because of the 5 changes in 
b 4 . This does not appear practical, inasmuch as it is a simple matter to change the 
load card at the CPC after the completion of the basic 10-group pass. In practice we 
have found it sufficient to number the groups up to the point where the first modification 
on the outermost loop is encountered. This modification may be thought of as the "major" 
change in the calculation. In our programming of the instruction deck we have assigned 
the set-number zero to Block 1, as at this point it is just a dummy. Later, at the CPC, 
this dummy will be replaced in turn by Block-sets 101, 102, 103, 104, and 105, which 
are shown programmed in the data programming. 

The group numbering which has just been described is the heart of the numbering 
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system. It is the group number which makes the cards for executing loops unique. It 
therefore appears that the numbering up to the group numbering, in particular the 
block numbering, is somewhat arbitrary. One limitation on block-numbering is 
obvious; namely, load blocks must bear a separate block-number, since in general 
they will also have set-numbers. Aside from this, it is perfectly legitimate to assign 
a single block number to a sequence of operations provided only that that sequence 
remain unbroken throughout the calculation. This is true even if the operations in the 
sequence are quite unrelated from the computational standpoint. Notice in our example 
that we could have abbreviated the block-numbering by assigning a single block number 
to Blocks 4 and 5, since these two always occur together. The same is true of Blocks 
2 and 3, and of Blocks 8 and 9. Whether to make this sort of condensation will depend 
on the particular problem. Often it is desirable to include more blocks than necessary 
to hold the sequence numbering to a minimum. Then if a programming error is 
discovered, the task of renumbering will not be so prodigious. 

When the group numbering is completed, the number of copies of each block-set 
required is entered on the left hand side of the CPC Deck Assembly Request. The 
operator then makes the required number by 80-80 reproducing and groups the copies 
according to the group numbering scheme. As' each group is formed, the group number 
is gang punched into the cards forming that group (into columns 9-10 in our systems). 
When the gang punching is complete all the cards are sorted on sequence number, set 
number, block number, and group number. The deck is then ready for the CPC. 

At the CPC 

The properly assembled CPC deck is capable of standing on its own two legs, so 
to speak. It can (and should) be sequenced every time it is run. It can also be 
sequenced by the tabulator while it is being run, if the equipment is available. It 
can be dropped and sorted back into order. We suggest that everyone making a deck 
in this manner drop it at least once to assure himself that this works. 

The important thing to note is that in any problem there need never be more than 
two types of cards at the CPC. One is the assembled instruction deck and the other 



may be load cards needed for making the "major change" in the calculation. By 
making these two decks on different type cards, e.g. , corner cut or colored, the 
problem of CPC operation is vastly simplified. 

There are, of course, many variations of this system. For example, the form 
of the system may be dictated by the number of card-numbering columns available 
on the instruction card. Another variation concerns further expansion of the 
instruction deck. If the expanded deck is too small to run conveniently, it may be 
reproduced 80-80 without further numbering. Of course the cards in the deck are no 
longer uniquely identified, and any sorting operation must be followed by a collating 
operation to remove the duplicates. 
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Finally, there are problems which have not been worked out fully. One has 
already been mentioned in connection with testing the program. The other concerns 
numbering summary cards which are to be used in further calculation. As time 
passes, it is hoped that these and other problems will be solved. 



FLOW CHART FOR 



CHART I 



Yto.bjJ-S {F(X ,bi)+4 [F(X, i b i )+ F(X,,bi)+ ...+ FfX 9l b i )]+i[F(X ii b i ) + ... + FtX e ,b i j]+F(X l() fi)} 




SET X k '0 
CLEAR STORAGE 
L,M,N. 



CALC. 
iF<X ,bj), 
ACCUM. IN L 



CALC. 



CALC. 
ftFlX k ,b,> 



T 





t 




ACCUM. 


a 

30 


F<X|<,bj ) 




IN M 



ACCUM. 
55 F(X k ,bj) 
IN N 



ACCUM. 
35 F(X| 0> bj ) 
IN L 



I 



CALC. 

y*(L) + 4(M)+2(N) 

PRINT y,b| 



T 



(3 




IS K-IO? 
YES NO 



IS K ODD? 
YES NO 



IS i-5 ? 
YES NO 



_T 



AEC-SE. RICHLAHO. WASH. 
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a Method for determining the roots of a secular determinant 
of order m, by means of the ibm model 602a calculating punch 

Andrew Acrivos and Neal R. Amundson 
University of Minnesota, Minneapolis 14, Minnesota 



I. Mathematical Apparatus 

In many applications of quantum mechanics and in a number of other problems it 
becomes highly desirable to be able to evaluate the roots of the secular equation D(X) = 
la tJ -bij XI = 0, i, j = l,2...m. 

To evaluate the m roots of the above secular equation, we will make use of the 
procedure shown in (1) and (4). This method consists of first reducing the secular 
determinant to a polynomial of the m* n degree in -A; the roots of this polynomial are then 
evaluated by means of any of the standard methods, i.e. , Newton's formula, Graeffe's 
formula, etc. 

Now, to transform the secular determinant into a polynomial in X one follows these 
steps: 

a. Let X - X = 0, 1, 2. . . m, thus obtaining m+1 different determin- 
ants; Xo is any convenient number. 

b. Evaluate these m+1 determinants. 

c. Then, by means of a simple procedure to be shown later, one 
uses these results to obtain the desired polynomial. 

Step b is by far the most troublesome, since the evaluation of an m* n order deter- 
minant is by no means an easy task, if m is larger than 3 or 4. It is here that the Model 
602A will be found useful. But first let us see how one is to go about evaluating a 
determinant. 

The simplest procedure consists in using the so-called "matrix inversion method", 
shown in detail in (2). We will give here a brief outline of this method. Consider for 
example the determinant which is obtained by letting X equal zero. 

a.u ais ais . . . aim 
A = a2i a22 a23 ... a2 m 



flml OjnS — in3 ... <*nun 
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Step 1: Divide each element of the first row by (-an); then write down the matrix 

-aia/an -ais/an . . . -aim /an 

10 

1 

♦ • • 

■ • • 

1 

This matrix has m rows, but m-1 columns 

Step 2: In the original determinant cross out completely the first row, thus obtaining 
the matrix 

&21 &22 a 2 3 • • • &2m 

a<n as2 as3 . • • a 3m 



axaJL am2 3m3 • * ■ «nun 

This matrix has now (m-1) rows, and m columns. 

Step 3: The two matrices obtained from the previous two steps are conformable, i. e. , 
one matrix has as many rows as the other has columns; therefore the matrices can be 
multiplied. So Step 3 consists in premultiplying the matrix obtained in Step 1 by that 
formed in Step 2. The result is the following: 

- a " a a i + a 2a - " a 2 i + a M ... 

an, a u 

- HZL a 3 i + a 8 2 - _fif a 3 i + ass ... 
an an an 



- JtOL a m i+ a^ - _*i! a^ + a m3 . . . - aim a m i + a m 

au an an 

It is not difficult to see that the above is again a determinant, but of the (m-l) th 
order. 

Step 4: Steps 1, 2 and 3 are now repeated over and over again until the original deter- 
minant is reduced down to one term, i.e. , a determinant of the first order. 

Step 5: The value of the original determinant is equal to the product of the left hand 
upper corner elements of the determinants obtained in the preceding steps, from order 
m down to order 1. Those are the same elements which are used, in each case, to 



aim 


aai + a 2m 


an 




ai m 


au + a 3m 
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divide the elements of the first row of each determinant, thus forming the matrices 
similar to the one obtained in Step 1. 

Let us now consider more carefully the nature of the arithmetic operations which 
are carried out in each one of the above steps. Step 1 clearly involves a simple division, 
while Step 3 requires an operation of the form 2 a„ • b n = P, where the summation 
includes only two terms. Step 4 is obviously a repetition of the operations carried out in 
Steps 1 and 3, while Step 5 consists of a series of multiplications. So one sees quite 
readily that the evaluation of a determinant of the m tn order can be done by means of 
only two arithmetic operations which are repeated a number of times. If we are then to 
use the Model 602A to evaluate said determinants we must be able to construct a wiring 
diagram which will enable the machine to perform said two operations, viz. , a division 
and an additive multiplication. 

Before proceeding with the wiring diagrams, however, let us say a few words 
concerning the degree of accuracy to which each of the two operations mentioned must be 
carried. One must bear in mind that due to the large number of operations involved, 
small inaccuracies introduced during the evaluation of the determinants might possibly 
be the cause of a rather serious error in the final result. We actually found by 
experience that, for determinants up to the 5th order, carrying one's calculations to the 
8th decimal place gave results that were at least within the desired degree of accuracy; 
i. e. , the computed roots were accurate to at least the 4th decimal place. Now it will be 
seen that the division, as set up on the Model 602A, carries the quotient up to the 8th 
decimal place only, and it is hoped that this will be sufficient for the accurate evaluation 
of a determinant of any order m. 

n 

n. The operation 2 aibi = P 

1:1 



We will try next to explain the wiring of diagram I which was set up to perform the 
above operation. Three types of cards are used throughout: 

1. The a card which carries an X in the first column. 

2. The b card which has an X in the 37th column. 

3. The P card which is blank, except for special cases 
mentioned below. 

The a card can have a 12-digit number punched in columns 5-16. The decimal 
point is between columns 7 and 8. The b card can have 15 digits punched in columns 
2-16 j the product P can have 14 digits punched in columns 3-16 and the decimal point is 
again between columns 7 and 8. If the numbers on a, b or P are negative, then an X 
must appear in column 17. 

The procedure used is quite easy to follow, inasmuch as it is standard. Since a 
has more than 8 digits, the operation has to be performed in 2 steps, as is shown in the 
diagram. 



-XtlXiRNATIONAl tUSINESS MACHINES COtfOS AT ION 

"ALCULATING PUNCH 







Panel Board I. Additive multiplication. 
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Now, since the answer card P has to be used in another operation, it must be 
punched as either an a or a b card. It the P card carries an X in column 59, then the 
answer card will be an a card, otherwise it will be a b card. Provision is also made to 
allow the machine to produce a certain number of copies of a given result. ..All copies, 
except for the last card, to be punched after each operation must carry an X in column 
77, to prevent the counters from resetting. 

The calculations are carried out to the 9th decimal place. 

s 

It should be noted that the sum 2 a t bi = P can contain an arbitrary number of 

terms. For more details, c.f. (3), pp. 17-19 and 52-53, where the simple multiplication 
is discussed quite extensively. 

HI. The Division (-b/a) 

The IBM Model 602A has been constructed to perform a simple division of the 
above form where the divisor a has only 8 digits. We are now going to propose a scheme 
by which the machine could be made to divide a number b, by a number a having as many 
as 12 digits. 

Let a = a' + c where a' has the first 8 digits of a while c has the last 4. Then we 
can write that 

u = iL( *_)= A (i-£_ + (.£_ + (-£-)"+...) 



a' + c a' 1+c/a' a' a' a' a' 

If now (c/a) < 10~* and if b/a' is not much larger than 1, then to the desired degree of 
accuracy, i.e. to the 8th decimal place, 

b b / a' - c v 
a" ss a 7( — 5 r_) 

In order, then, to evaluate the quotient b/a, if a has more than 8 digits, we would have 
to perform 4 successive operations, viz. , a simple division b/a', a subtraction of two 
numbers a' - c, another simple division (a' - c)/a', and finally a multiplication. As 
would be expected^ the wiring diagram for such an operation is quite complex, and it is 
for this reason that we are going to analyze it quite thoroughly in the following pages. 

Planning Charts and Control Wiring for the Division - b/a 

(Diagram II) 

a) a' is read into Storage Entry 1R while c is read into Storage Entry 2R If a is 
negative, i.e. if an X is punched in column 17, then pilot selector 3 is picked 
up. The a card is next skipped out and the machine is ready to read the b card. 
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INTERNATIONAL BUSINESS MACHINES CO«P£ 

CALCULATING/FtJNCH 
T Y fe aoa A -cgyrnoL panbl 




Panel Board II. Division 
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b) b is read into counters 1, 2 and 3 which are connected. If b is negative then pilot 
selector 4 is energized. 

Program Steps; 

Programs 1, 2, 3 and 4 perform the division b/a' carried out to the ninth decimal 
place. The operation is standard; c. f. (3), pp. 54-55. 

Program 5 is used to skip programs 6,7,8 and 9. 

Program 10. The quotient b/a' is read into storage entries 6R and 6L. 

Programs 11 and 12. c is subtracted from a' while at the same time pilot 
selector 6 is picked up. Counters 4-6 are reset to 5. 

Programs 1', 2', 3' and 4' perform another division, viz. (a' - c)/a\ The 
quotient is obviously less than unity and will contain 8 digits. The last digit has been 
reset indirectly during Program 11. 

Program 6. (a' - c)/a' is read into 1R while the counters are reset. 

Program 7. b/a' is multiplied by (a' - c)/a\ At the same time all the counters 
become connected by means of co-selector 8. 

Program 8. The final answer is read into storage unit 7 and the machine is 
ready to read the next card. 

c) The answer card is blank except maybe, as will be seen later, for an X punched 
in column 50. 

Program Steps: 

All program steps up to 9 are skipped. 

Program 9. The answer is punched out of storage punch exit 7; the counters are 
reset. Since the quotient is equal to (- b/a), the answer card must have an X punched 
in column 17 unless b and a have different signs. It is easily seen that if either pilot 
selector 3 or 4 is energized, then pilot selector 7 will be energized before punching and 
no X will be punched in column 17 of the answer card. One can, of course, reproduce 
any number of "copies" of the answer. It is very important, however, to note that the 
last answer card after each operation must have an X punched in column 50, otherwise 
pilot selectors 3, 4 and 7 will not drop out. 

In this particular case, an X is punched during the punching in the first column of 
the answer card; therefore the answer card becomes an a card. Of course, one could 
easily change the diagram so that the answer card would become a b card. 

Important Note: Let us for a moment consider one of the four operations which had to be 
carried out during the division (-b/a), i.e. (a' - c')/a'. The quotient is a fraction which 
comprises 8 digits. Those 8 digits are then used to multiply b/a' and thus produce the 
final answer. Suppose now, however, that c'-*- 0; in that case obviously (a' - c*)/a '-*•!, 
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and the 8 digits mentioned before turn out to be all equal to zero. If then one uses the 
wiring diagram unmodified, one obtains an answer equal to zero instead of b/a. It is 
desired to show here that the diagram can be altered very easily to perform this simple 
division. 

So if c = 0, the wiring diagram should be changed as follows: 

1. Disconnect the "Read" from program 8 and connect it to program 5. 

2. Read out of 6 and into Storage Unit 6 at program 5, instead of pro- 
gram 10. 

3. During program 9, read out of 6 and into Punch 6 instead of reading 
out of 7 and into Punch 7. 

4. Connect column 16 in the Punching section to the first hub of Storage 
Punch Exit 6, and the other hubs of Storage Punch Exit 6 to Punching 
Columns 5-15. Connect column 4 to the punch emitter. 

The machine then will give the answer of (-b/a) to the ninth decimal place. 

Limitations: The quotient (- b/a) will be correct to the 8th decimal figure if: 

1. b/a' is not much larger than one. 

2. If (c/a') < NT 4 . 

Also a can have at most 12 significant digits (9 decimal places) and at least 9 
significant digits. The wiring diagram must then be altered if it is to serve other more 
specific cases. 

IV. How to Evaluate the Roots of an m Order Secular 
Determinant by Means of the IBM Model 602A 

Step 1: Let \ - X = 0, 1, 2, . . . m* and set up the m+1 determinants thus obtained. 
The term an must be an a card, while all the other elements of each determinant must 
be b cards. 

Step 2: Obtain, by using Panel Board n, the m+1 matrices analogous to the one 
shown in Step 1 of Part I, performing the divisions continuously (a certain number of 
copies of each answer must also be produced by the machine). 

Step 3: Use is now made of Board L Follow the procedure explained in Step 3 of 

Section I. We thus obtain m+1 determinants of the (m-l)th order. It must be noted that 

the term r - ai2 a^ + a 22 i (left hand corner) must be punched as an a card while all the 
a u 

other elements of the determinant must be b cards. In all cases, the additive multipli- 
cation contains two terms only, i. e. , it is of the form ef + g=p; this can also be 
written as e f + 1 g = p, which, of course, is the operation Board I is supposed to perform. 
In the same way, i. e. , by repeating Steps 2 and 3 over and over again, one can 

* For (X - Xo) = h, 2h,...mh where h* 1, c.f. (4). 
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evaluate the m+1 determinants, all at the same time; this is explained in Section I. 

Step 4: With the values of those m + 1 determinants at hand, one can easily reduce 
the secular determinant to a polynomial of the m degree on X. 1 (Consult (1) and (4).) 
This polynomial is of the form 

f (X) » X m + aX- 1 + bX— 2 + . . . + yX+w 

and 

f '(X) = mX*- 1 + a'X"- 2 + . . . + y 

where f ' (X) is the derivative of f (X) with respect to X. The Panel Boards I and II can 
now be used to compute the roots of f (X) by means of Newton's method, if all the roots 
are real. The following procedure is recommended: 

Stepl: Use Board I and evaluate f (X) and f ' (X) for an assumed value of X. 

Step 2: By means of Board II, compute the quotient ■- f (X)/f'(X). 

Step 3: A better approximation for one of the roots will be given by 

X' = X + (-*<*) ). 
f'(X) 

Repeat Step 1, 2 and 3 until all the roots are determined to the degree of accuracy 
desired. 

V. Conclusions 



It must be borne in mind that the method, as outlined above, is quite simple in 
that at no time during all the various operations one is required to read any of the 
intermediate answers obtained. Only the final answer, i.e. , the X's, should be read out 
of the cards. This is a distinct advantage, as it eliminates human errors and saves time. 
Moreover, aside from the punching machine, the only other equipment used is the Model 
602A Calculating Punch, which performs all the operations required. Finally, the 
evaluation of the determinants can be checked by a simple procedure mentioned in (1) 
and (4). 
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NOTES ON THE 602-A AND 604 
(Supplement to Technical Newsletter No. 3) 

Eric V. Hankam 

Watson Scientific Computing Laboratory 
International Business Machines Corporation 



1. The Balance (punch) conversion unit can only be activated by a or 9. Any other 
digit wired into the B. T. hub will punch in the high order position in place of the 
actual first digit. 

2. To p.u. a selector at the end of P12 (when more than 12 steps are required) and if 
MPL or DIV takes place on P i2 , wire: 




Pi P12 

(Supplement to Technical Newsletter No. 3, p. 56) 

3. a) Can multiply on two successive steps by X2 or X5; also ordinary multiplication 

can be followed by X2 or X5. 

b) Cannot MPL or DIV on Read cycle. 

c) Impulsing MPL and TENTH will not shift product one place to the right. 

4. Emitter (but not Digit Selector converted into emitter) can be used as a constant 
multiplicand. 

(Supplement to Technical Newsletter No r 3, p. 54) 

5. Any wires from the Emitter (not Digit selector converted into emitter) to a counter 
entry should be isolated during MPL or DIV steps (even if that counter is not involved 
in the MPL or DIV operation, and even if that counter is not impulsed to RI during 
the MPL or DIV step). 

6. Whenever two (or more) counter exit positions (of the same counter or counter 
group) are split wired together, that counter should not be reset on the same step on 
which it is RO, and no other counter should be RI on the same step. An internal 
emitter makes common all like digits (on a given program step), and therefore also 
unlike digits if any two of them are split wired. 
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7. It is possible to RO one counter of a group of counters that are coupled together, 
and, on the same step RI another counter of the same group. In particular, the 
number in one counter of the coupled group can be transferred to another counter of 
the coupled group. 

8. Whenever a counter is simultaneously impulsed to RI+ and RI-, it will not read in 
at all. 

9. Division with a complement dividend can be done as described in Technical News- 
letter No. 3, p. 59, but if the division requires 7 shifts, the first division cycle 
must be wired differently as outlined in the diagram below.* With this modification, 
any division problem with a complement dividend can be handled regardless of the 
number of shifts required. 

X-p.u. 



I-p. u. 



[1,2,3]RI- 




[7,8]RI+ 
(quotient) 



hot X from 
Digit Emitter 



Program exits 



^ 



lR)RO 



* This method was suggested by Mr. S. Kelly, IBM Customer Engineer. 
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16 ± 16 ADDITION ON THE 604 
(Supplement and correction to Technical Newsletter No. 3, p. 50) 



a) On machines equipped with Zero Test, a + bw-o can be zero tested directly on P 6 , 
and several of the subsequent program steps can be eliminated. The a Adjust 
operations on P 19 and Ps* will be suppressed if a + bn-e = 0. 

b) The £ Adjust operations on P 19 and P^ should also be suppressed if a + bis -u* 0. 
For that purpose, program steps should be inserted to zero test a + bia-». The 
number in EC has to be "boosted" by ±5 to preserve its sign only if a + bia-o is zero, 
but this method would fail in case a + bn-e= 9995. 
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15 X 15 MULTIPLICATION ON THE 604 

Eric V. Hankam 

Watson Scientific Computing Laboratory 
International Business Machines Corporation 



I h 

a 

I 
Fl,3 



b 
F2 



c 

i 

MQ 



A 
Gl*3 



B 



C 



H 



G2 free 



PI 

P2 

P3 

P4 

P5 

P6 

P7 

P8 

P9 

P10 

Pll 

P12 

PI 3 

P14 

PI 5 

P16 

P17 

PI 8 

P19 

P20 

P21 

P22 

P23 

Punch 



cB 

c — 
b - 
bC 

c — 
A — 
M 4 - 
M 4 - 
cA 

C - 
aC 

B - 
bB 

round 
M 3 — 
M 3 - 
aB 

A — 
bA 

M 2 — 
M 2 — 
M 2 — 
aA 



G2 
MQ 

MQ 

G2 

Gl,3 

EC 

MQ 
MQ 



Gl,3 

EC 

MQ 

G2 

Gl,3 

EC 



EC 

4 

15 positions 



G2 

5 



F4 RO ; MPL+ =► 

MQ RO ; G2 RI 

F2 RO ; MQ RI 

G4 RO ; MPL+ ==► 

G2 RO ; MQ RI 

Gl,3RO ;G2RI 

ECRR;G1,3 RI; Ro 6th 

Gl,3 RO ; ECRI+ 

G2 RO ; MPL+ => 

G4 RO ; MQ RI 

Fl,3RO;MPL+ ==>. 

F4 RO ; MQ RI 

F2 RO ; MPL+ =► 

1/2 Adj ; Ri 5th 

ECRR ;G1,3 RI ; Ro 6th 

Gl,3RO ; ECRI+ 

Fl,3RO;MPL+ => 

G2 RO ; MQ RI 

F2 RO ; MPL+ ==*• 

EC RO ; G2 RI 

EC RR;G1,3RI; Ro 6th 

Gl,3 RO ; EC RI + 

Fl,3RO;MPL+ >. 



cB e EC 



cB + bC = M 4 € EC 



(highest 6 places) 

M 4 + cA e EC 

M 4 + cA + aC e EC 

M 4 + cA + aC + bB = M 3 e EC 

(highest 6 places) 

M 3 + aB e EC 

M 3 + aB + bA = M 2 e EC 
(lowest 5 places) 
(highest 6 places) 

M 2 + aA = Mi e EC 



NOTE: This method differs from a previous article (Technical Newsletter No. 3) in 
that it includes the partial products cB and bC , thus reducing the round-off error. 
The previous method has the advantage of requiring fewer electronic cycles. 
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16 x 16 MULTIPLICATION ON 604 

Eric V. Hankam 

Watson Scientific Computing Laboratory 
International Business Machines Corporation 



3 5 

I h 



b c 
k G2 



k A * 



3 

-t 1 



d 



3 5 5 3 

i 1 1 1 

A B CD 

F3 F4 MQ Gl 



\ 



G4 free 
(8 - 6) assignment 



PI 


: C-*G4 




P2 


: C-»EC 




P3 


: C-»MQ 




P4 


: c C 




P5 


: c C-*MQ 




P6 


: EMEC 




P7 


: c C->G1 




P8 


: D-»MQ 




P9 


: c C^EC 




P10 


a-*Gl 




Pll 


(a + B) D 




PI 2 


A-K51 




PI 3 


d-»MQ 




PI 4 


(A + B) d 




PI 5 


I M x 


MQ : 


P16 


J """-> 


Gl : 


P17 


1 






\Mx »• 


EC 


P18 : 


J 




P19 : 


c-*MQ 




P20 : 


Be 




P21 : 


C f 


MQ : 


P22 : 


bC 





P23 



round 



MQ RO ; G4 RI 

MQ RO ; EC RI+ 

EC RR ; Ro 3rd; MQ RI ■ 

G2 RO ; MPL+ 

EC RR ; Ro 6th; MQ RI = 

GIRO ; EC RI+ 

MQ RO ; Gl RI; Ri 6th 

EC RR ; Ro 6th; MQ RI '■ 

Gl RO ; EC RI+ 

Fl RO ; Gl RI; Ri 6th 

Gl RO ; F2 RO; MPL+ = 

F3 RO ; Gl RI; Ri 6th 

G3 RO ; MQ RI 

GIRO ; F4RO; MPL+ 

EC RO ; MQ RI; Ro 5th 

EC RR ; Gl RI; Ro 5th 
MQ RO ; EC RI+; Ri 2nd 

Gl RO ; EC RI+; Ri 2nd 

G2 RO ; MQ RI 

F4 RO ; MPL+ = > 

G4 RO; MQ RI 

F2 RO ; MPL+ 

1/2 Adj. ; Ri 5th 



• C € MQ ooxix dropping lowest 2 

C C € EC xxxlcxxxx 

highest 3 pos. of c C in MQ 

xxx ooooo (due to 8 - 6 assignment) 



•D e MQ ooxxx 



XXX ooooo 



xxx ooooo 



cC 



c C + (a + b) D e EC * x x x\ x x*x*x7x * 



xxx ooooo 



cC + (a + b) D + (A + B) d = Mx e EC ^i,,,,,,, 

v ^ — *\ 



000 xxxxx 



XXX OOOOO 

o oxxxxx o 



Mi + Be € EC 



^— -C a r r y 
if o xxxxxxxx 

X X XX X XX XX X X 



Mi + Be + bC = M 2 e EC xxxxxxixxxxx 

t 
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P24 

P25 
P26 

P27 
P28 
P29 
P30 
P31 
P32 
P33 

P34 
P35 
P36 
P37 
P38 
P39 

Punch 



} 



) 



M 2 

M 2 

Ca 

c 

cA 

B 

bB 

M 3 

M 3 

M s 
aB 
A 
a 



G4 

Gl 
EC 



-+• MQ 
-* MQ 



-> G4 

„ G2 

"* Gl 

-*- EC 



-> MQ 



Gl 



(a + b) A 




EC RO ; Ro 6th ; G4 RI 

EC RR j Ro 6th ; Gl RI 
Gl RO ; G4 RO ; EC RI+ 

Fl RO ; MPL+ 
G2 RO ; MQ RI 
F3 RO ; MPL+ 
F4 RO ; MQ RI 
F2 RO ; MPL+ 
EC RO ; G4 RI 
EC RO ; Ro 6th ; G2 RI 

EC RR ; Ro 6th ; Gl RI 
Gl RO ; G2 RO ; EC RI+ 
Fl RO ; MPL+ => 
F3 RO ; MQ RI 
Fl RO ; Ri 6th ; Gl RI 
Gl RO ; F2 RO ; MPL+ = 



=> 



o ozo oooo 



xxxxxx 



+ Ca € EC xxxxxxxxx' 



.M 2 



M : 



M 2 +Ca + cAeEC ,„„„„. 



M a + Ca + cA + bB = M 3 e EC xxxxx"""" 1 
(lowest 5 digits) to be punched 



OO OX X XXX 



OXXXXX A 

xooooo J 



highest 6 digits 



xxxxxx 



(highest 6 digits) 

M 3 + aB=M4€EC ,„,„„ 



►M 4 + (a + b)A € EC 



xxxxxxxx <- 
xxxxxxxxxxx. 



.M4 



11 pos. 



16 positions 



NOTE: Read in and calculate absolute values; determine final sign through selectors. 
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DIFFERENCING ON THE CPC - MODEL II 

Stanley Poley 

Department of Civil Engineering 

Columbia University 

New York, New York 



A description is given herein of a set of control panels for the CPC - Model n * 
by means of which high order differences may be easily and rapidly computed. 
Principal features of this set-up are: 

1) The order of the highest difference desired is set by plugging one 
wire. 

2) No time is wasted in computing differences of order higher than 
called for in (1). 

3) No blank cards or instruction cards are required; only function 
cards are run through the 418. 

The control panels described in Part I will accommodate differences up to tenth 
order for ten digit numbers. In Part II, several auxiliary features are given which 
increase the capacity to twentieth order differences of thirteen digit numbers. 

These control panels were wired and tested at the Watson Scientific Computing 
Laboratory at Columbia University. 

Parti 

The flow of data is shown in Figure 1. Note that one 418 card cycle and n 418 
program cycles are required to obtain an n* n order difference. Thus the speed of 

of operation is =- cards per minute. 

n + 1 



!• Equipped with a special program unit. 
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Cycle 


Card 


6A 
6B 


941 Storage Units 


F2,4 


G2,4 


EC 


List 
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12 


13 
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_\7 B 
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V}° 




P10 

.... 














T7 9 £ 






RO— 


»V}-° 




Vi* 







Figure 1 : Flow of Data 



Wiring for the 418 is shown in Figure 2, while Figures 3 and 4 give, respectively, 
605 program pick-up wiring and 605 programming. 

PROG. SOURCE 



T 0- 



N ©- 



Cjp 



CAL-SEL 1 



PROG. PU 



->Q-A=2 — O 



->Q 4 - 6 O 



Figure 3: 605 Program Pick-up Wiring 



Prog. 


Step 


Operation 


1 




G4 RO, EC RI+ 


2 




G2 RO, EC RI+, Ri 6th 


4 




EC RO, G4 RI 


5 




EC RO, G2 RI, Ro 6th 


7 




F4 RO, EC RI- 


8 




F2 RO, EC RI- , Ri 6th 



Figure 4: 605 Programming 



All differences are stored in the 941 and are read in or out by emitting channel 
codes through successively transferred levels of a field selector, level corresponding 
to PI, level 1 to P2, . . . , level 9 to P10. However, the function y^ in Figure 1 is 
stored in 418 counters 6A, 6B in order to avoid 941 RI-RO timing difficulties which 
would be encountered when taking first differences. 

As indicated in Figure 2, one wire must be plugged to level (n - 1) of the right- 
hand column of the field selector to obtain n th order differences. This wire impulses 
MI-STOP and transfers Co-Sel 55. 2 Since the field selector is picked-up through the 
normal side of Co-Sel 55, it will not be picked-up during Pn and will thus be normal 
on the following card cycle. Channel A control is wired through the normal side of 
Co-Sel 55 to prevent a 941 RO during card cycles which would conflict with the RR of 
counters 6A, 6B. 

2 Notation in this paper is that employed by Eric V. Hankam in his paper, "Notes on 
the IBM Type 604 Electronic Calculating Punch and Type 602A Calculating Punch", 
IBM Technical Newsletter No. 3. 
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The card column containing the minus punch for y is wired from 2nd reading to 
the common of a digit selector. Pilot selector 70 should be picked-up from the digit 
selector hub corresponding to the minus punch. 

The differences list from EC through five gating selectors whose pick-up hubs are 
denoted by A, B, C, D, E in Figure 2. These pick-up hubs should be wired from 
correspondingly labeled hubs of the field selector. When gating selector E is trans- 
ferred, Pil-Sel 48 will transfer on the following cycle, suppressing space suppress for 
that cycle. Thus the differences list in two rows, five across each row. The function 
y lists to the left of the first row. 

Part II 

A) Extension for Higher Order Differences: 

By wiring a second field selector in series with the first, as shown in Figure 5 
(ordinary selectors may be wired to operate like a field selector if only one field 
selector is available), up to twentieth order differences (of ten digit numbers) may be 
obtained in a single run. 

B) Thirteen Digit Numbers: 

By wiring 418 counters to operate in conjunction with the 941 storage units, we 
may accommodate up to thirteen digit numbers if it is assumed that differences of order 
greater than eleven will not exceed ten digits. The additional wiring required is shown 
in Figure 6 (for up to twentieth order differences). Moreover, Fl and Gl should be 
given (8-6) assignments and the 605 controls of Fl, F2 and Gl, G2 coupled, i. e. , in 
Figure 4, F2RO becomes Fl,2RO, etc. Counter wiring associated with the second field 
selector may be omitted if only tenth order differences or less are required. 

C) Designation of "n" by External Switches: 

In some installations, it may be desirable to select n , the order of the highest 
difference, by using external switches rather than by plugging wires on the 418 control 
panel. Figure 7 shows a Christmas tree circuit wherein n is set in binary form using 
AB switches (set-up change switches may be substituted for AB switches, if necessary). 
Setting of the switches for n = 19, for example, is shown in Figure 8. Note that the AC 
switch position corresponds to 1 while BC corresponds to . 
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Emitter 



Pil-Sel , T n Pil-Sel 

56 ° N ° 54 



Digit Pulses from Emitter 
Field Selection 



-Select "n' 




MI-STO 



• 10 Common 



A 
Channel 
Control 



FSPU 



Figure 5: Field Selector Wiring for up to 20th Order Differences 
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Figure 7: Wiring for External "n" Designation 
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Figure 8: Switches Set for 19th Order Differences 
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CALCULATION OF HIGH ORDER DIFFERENCES 
ON THE 607 ELECTRONIC CALCULATOR 

Eric V. Hankam 

Watson Scientific Computing Laboratory 

International Business Machines Corporation 



On a high speed calculating machine equipped with Program Repeat, such as the 
607, where the skipping of program steps does not involve any appreciable loss of time, 
certain types of problems characterized by a repetitive pattern can be set up by repeat- 
ing a "loop" of relatively few program steps. This set of program steps can be con- 
sidered as a "subroutine" which can be incorporated on a control panel on which many 
other operations are wired. The evaluation of high order differences has been chosen 
to illustrate this principle. 

The program below describes the computation of 8th order differences of 8-digit 
numbers. This can easily be extended to a) higher order differences, and b) 13-digit 
numbers. If Program Repeat Selectors were available, a selector network could be 
set up to count the number of program repeats, and terminate the calculation after the 
required number of repetitions. Such a selector network would further decrease the 
number of program steps in the loop. In the example following, steps 13-15 perform 
the counting by successive subtractions of l's from 7 (one less than the highest order 
difference to be computed). Group Suppressions might provide an alternate method to 
keep track of successive program repeats. 



Read: 



fi 
7 



607 Program for 8th Differences* 



-> Fl,2 

-*- F3 



PI 
P2 



Fl,2RO 
ECRO 



; EC RI+ 
; G17,18RI 



P3 
P4 



Gl, 2 RO ; EC RI- ; SKIP OFF (develops Af , k = 1, 2, . . . 8) 

G3,4RO ; G1.2RI 



,-th 



♦Description is for i in card. Notation used is similar to that employed in Technical 
Newsletter No. 3. 
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P5 


G5,6RO , 


G3,4RI 


P6 


G7,8RO 


G5,6RI 


P7 


G9, 10 RO 


G7,8RI 


P8 


Gil, 12 RO 


, G9, 10 RI 


P9 


G13,14RO 


, G11,12RI 


P10 


G15,16RO 


G13, 14 RI 


Pll 


G17, 18 RO 


; G15,16RI 


P12 


EC RR 


; G17,18RI 


P13 


F3RO 


; EC RI+ 


P14 


: Em 1 


EC RI- ; B. T. 


PI 5 


EC RR 


; F3 RI 


P16 


: G17, 18 RO 


; EC RI+ ; Prog. Repeat ; Supp 


PI 7 


SKIP ON 





Punch: At the end of the calculation the differences will appear in and can 
be punched from the following units: Aj in G3, 4; A? in G5, 6;. . . 
A? in G17, 18. The differences remain in these storage units for 
the (i + l)st card. 



f3 



Schematic Diagram of Program Loop 
fi 



^ EC (Electronic Counter) 



T>Z, 7>/z 




Note that the (i - l) st 8th difference is first replaced by the 1 th function on P2 , then 
the i" 1 first difference is obtained on P3 , and then each difference is successively 
shifted one storage unit to the left, etc. 

Remarks on the Program 

1 • a) f i need not be read from the card; it can be the result of a preliminary compu- 
tation, and must be entered into EC before the differences can be computed. The 
7 can also be emitted into F3 on any program step prior to the computation of the 
differences, 
b) If fi is read from the card it could be entered directly into G17, 18 instead of 
Fl,2; in that case A 8 may be transferred to EC to be punched. 

2. The result of the B. T. on Pj^ can be used to control program steps which are not 
part of the differencing loop, e. g. SKIP ON on Pj 7 can be suppressed on - if 
other calculations are to be programmed after the differences have been obtained. 
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RAPID DIFFERENCING ON IBM ACCOUNTING MACHINES 

E. Brown 

Physics Dept. , Cornell University 

Ithaca, New York 

Introduction 

A standard operation on IBM accounting machines is the calculation of differences 
of a function. The methods currently in use require program cycles on the accounting 
machine which may add considerably to the calculation time. The method discussed 
here is one in which program cycles are not required. 

Method 



The basis of the present method is in the use of a group of three counters to obtain 
the first differences of the function read into them. These counters are instructed in 
accordance with the following scheme: 



Machine Cycle 

# 


Counter 
Ai 


Counter 
A a 


Counter 
A 3 


3n 


RI+ 


RI- 


Total 


3n+l 


Total 


RI+ 


RI- 


3n+2 


RI- 


Total 


RI+ 



It can readily be seen that, if successive values of a function are read in on each 
cycle, the differences will be successively developed in each one of the three counters. 
The difference A^f^-f^ will be read out of one of the counters on the same cycle 
that f n+a is read into the other two. If the output of this group of counters is read into 
a similar group of counters, then the output of this latter group is the second difference 
of the function. This procedure can, in theory, be used to obtain any order difference. 
The time required is independent of the order. The restriction on the order of the 
difference is imposed by the number of counters in the machine. 

The panel described in this paper obtains third order differences of seven digit 
numbers and makes use of nine eight-position counters. (A combination of a coupled 
6 and 2 position counter will be referred to as an 8 position counter). These nine 
counters will be referred to by the following notation: 
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Ai, A 2 , A 3 , 
Bi, B2, B 3 , 
Ci, C2, C3, 



-> counters necessary for first differences 
-> counters necessary for second differences 
-^ counters necessary for third differences 



The counters are controlled by means of two pilot selectors. The selector wiring 
is illustrated below: 




Card 
Count 

Assuming both these selectors to be normal to start with we obtain the following 
pattern: 



Machine Cycle # 


Selector # 1 




Selector # 2 


3n 


N 


N 


3n+l 


N 


T 


3n+2 


T 


N 



The counter controls are impulsed through the other position of the pilot selectors 
as shown below: 

Sel. #1 Sel. #2 

N ^ N 




RI+ 



Total 



RI- 



(|) EI - (1)™*) 

(th (lW 



For simplicity in wiring, the exits of the counters in each group are connected. 
Because of the fact that a hot pulse reaching a counter exit will prevent it from 
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resetting properly when it is impulsed to total it is necessary to impulse counter exit 
suppression for each counter when it is reading in (plus or minus). 

The method illustrated here yields true figures for the differences. For this 
purpose the NB exits of the counters are wired to the corresponding NB control entries. 
Thus, when a counter which contains a negative number in 9 complements form is 
totalled a conversion cycle takes place which transforms this number to a true figure. 
These conversion cycles add time to the calculation. (It may be possible to remedy 
this defect by using the usual procedure for subtracting on non-net balance machines. 
This procedure however requires twice as many counters as are required here. ) How- 
ever, if one is willing to print complements, one may obtain the maximum speed of one 
card cycle per set of differences. The wires which are represented by dashed lines 
may be removed for this purpose. For greater flexfcility it is preferable to wire them 
through co-selectors which may be picked up from a set-up change hub. 

The purpose of the field selector and associated selectors is to prevent printing of 
meaningless differences at the beginning of a run. For example the first difference 
column has no entry until f s is listed. The printing takes place as follows: 



Function 


First Difference 


Second Difference 


Third Difference 


fi 








u 








fs 


fa-f!= AJ 






u 


f 3-fa= A* 






u 


I *-Uf A J 


Aj-AJ=A? 




u 


f 5 -fr*A$ 


AJ-A*=A§ 




f 7 


fs-fiP M 


AJ-4=A| 


Ai-Aj=Af 


fs 


h-i e= Ai 


A*-Ai=A 4 * 


As*-A§=A! 


• 


• 


• 
• 


• 



The card deck should be followed by three blank cards in order to obtain the 
last significant third difference. This will result in the printing of a meaningless 
second difference for the last card and a meaningless first difference for the preceding 
card. No provision has been made to suppress these on the present panel. If one 
desires to suppress these differences from printing it may be achieved by usingx-punches 
in the trailer cards to pick up the co-selectors through which these print. 

The panel presented here was wired and tested at Watson Scientific Computing 
Laboratory, International Business Machines Corporation. 
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Notes 

1) It is possible for the field selector to be transferred at the beginning of a run. This 
is the case if a deck has just been run since the LC pulse initiates the transfer of the 
field-selector. For this reason the power switch should be turned off and on prior to 
running a deck. 

2) It is not necessary to clear any counters before starting. 

3) The function cards should not be preceded by a blank card. 
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ERRATA 

In Technical Newsletter No. 5 in the article "Boolean Notation for Selector Net- 
works" by Eric V. Hankam, there is an error in the diagram (d) at the top of page 27. 
The correct diagram should be: 

T 

N / „ A 1 v A J 

C 




Technical Newsletter No. 5 - A Model H Four-Address Floating-Decimal Coding 
System - Dura W. Sweeney. 

Mr. Richard C. Lesser of Cornell University has discovered the following three 
inaccuracies which appear in the calculate portion of the CPC board described in this 
article. 

1. The computation of e where 

B » -n /n 10 

where n = 1, 2, 3, 4, 10, will cause the 
machine to compute as if 

B = n Jfn 10 

This occurs since division by Xn 10 now gives a minus zero remainder, and on 
a later program this sign is lost as the remainder is added into the electronic 
counter. Since n is later half-added into the counter under control of the sign 
of the remainder, this causes the answer to appear as 10 rather than 10 . 

2. When the number 9. 9999995 is added to zero (or, in fact, if 9. 9999995 is the 
result of any operation), the machine, by virtue of double round-off, gives 
10. 000000 as the answer. 

3. Although the manual (LA 1475) states that multiplication can be performed if 
numbers are not in true floating decimal form, a restriction regarding 
significant figures should be noted. 

Multiplying 00741000 x 00741000 

0. 0741 000 54 741. 

0. 0741 opq 54 741. 

0000549081000000 549UBTT 

the machine rounds-off after the 8, and the 1 is dropped, giving, after shifting 
a result of 

5.4908000 55 



78 



• half-add 


o P 


o ctr + 


• #15 


o (63) 


#26 



Mr. Sweeney has made the following remarks and modifications in connection with 
these. 

1. I have modified the 605 control panel as follows to eliminate the difficulty of 
obtaining 

e -n/nl0 

1. Panel C step 49 change "half-add" to "ctr + (63)". 

2. Panel C step 51 change emit 2; divide; — ; #44 to RO; GS3; — ; #15. 

3. Panel C step 52 change emit 1; ctr -; BT; reset; — ; #44 to emit 2; divide; 

— ;#34. 

4. Panel C step 53 change RO; GS3; — ; # 15 to emit 1; ctr -; BT; reset; 

II A 1 

— ; it ot. 

5. Add to explanations of selected operations: 
(63), CS 11, half-add (S), ctr + other. 

6. Change CS11 from 

To © P 

N • #44 « SUP to 

C o SUP • #26 

2. The difficulty of rounding off . 99999999 >x >. 99999995 to 1. 0000000 is a 
compromise. Only if a result is in this range is it so rounded. If this rounding 
were eliminated, all other operations which "overflow" would be incorrectly 
rounded. Since these few cases are an extremely small percentage of the 
possible results, it does not seem worthwhile to change the control panels. 

3. The third case mentioned is inherent in the design of the panel. The case of 
inaccuracies in multiplying non-shifted numbers can also occur in division or 
addition. For example: 

00000234 50 + 01791342 50 = 13062832 46 

actually gives: 13063000 46 
00000234 50 + 17913421 45 = 41313421 45 

actually gives: 41300000 45 

It should also be noted that in the operation (A + B)/ X, the result of the addition 

is not shifted before the division takes place. Specifically, for A = 12345876, 

B = 12345731, X = 7, the result of the above operation is 207 rather than 20714286. 

Technical Newsletter No. 6. 13 Digit Floating Decimal Model II CPC - Dura W. 
Sweeney. 

Dr. F. M. Verzuh of M. I. T. brought the following inaccuracy in this description to 

our attention. 

In using this panel to perform the following division: 

Q = A = 4.472135955000 
B 8. 944271910000 
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An incorrect quotient of 0. 4999899999954 is obtained instead of the correct value of 
0.5. 

Mr. Sweeney has suggested the following changes in the Panel A wiring to eliminate 
this error. 





PANEL A 




Prog. 




Exit 3 


17 




In 5 


18 




In 5 


19 




In 4 



27 mi 

28 Out 2 



30 In 5 

In the same article the 605 Planning Chart under Panel B, page 39, indicates 
that Repeat Selectors #3, 4, and 10 are to be picked up by exit 1 of Program Step 59. 
Experience has proved that multiple wiring of Repeat Selectors has been extremely 
troublesome. For this reason the following changes to Panel B are suggested. 

PANEL B 

Prog- Exit 1 Exit 3 

47 Prog. Rpt. , Rpt. Del PU RS # 10 PU 



59 RS#3, 4 GS1&2RI 



Technical Newsletter No. 5, p. 10 line 10 read 1570796350 for 1570596350. 
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